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Abstract 


We continue to study the holographic p-wave superconductor model in the 
Einstein-Maxwell-complex vector field theory with a non-minimal coupling be- 
tween the complex vector field and the Maxwell field. In this paper we work 
in the AdS soliton background which describes a conformal field theory in the 
confined phase and focus on the probe approximation. We find that an applied 
magnetic field can lead to the condensate of the vector field and the AdS soliton 
instability. As a result, a vortex lattice structure forms in the spatial directions 
perpendicular to the applied magnetic field. As a comparison, we also discuss 
the vector condensate in the Einstein-SU(2) Yang-Mills theory and find that 
in the setup of the present paper, the Einstein-Maxwell-complex vector field 
model is a generalization of the SU(2) model in the sense that the vector field 
has a general mass and gyromagnetic ratio. 
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1 Introduction 


The gauge/gravity duality [1, 2, 3] allows us to map a strongly interacting quantum field 
theory to a weakly coupled gravity theory living in a higher dimensional spacetime. The 
applications of the gauge/gravity duality have been extended to various fields in recent 
years, especially in condensed matter physics [4, 5, 6]. 

A well known example is the holographic superconductor model [7, 8]. By introducing 
a charged scalar field coupled to a (3 + 1) dimensional Einstein-Maxwell theory with a 
negative cosmological constant, this model admits a second order phase transition. At 
high temperature, the scalar field vanishes and one obtains the standard AdS Reissner- 
Nordstrom black hole. However, when one lowers the temperature, the AdS Reissner- 
Nordström black hole becomes unstable and is replaced by a charged black hole with 
a non-trivial scalar “hair”. The higher temperature AdS Reissner-Nordstrém black hole 
describes a conductor phase, while the low temperature black hole solution has the expected 
behavior for a superconducting phase. Since the condensed field is a scalar field dual to a 
scalar operator in the field theory side, it is therefore a holographic s-wave superconductor 
model. Holographic p-wave models were realized by the condensation of a vector like 
field [9, 10, 11, 12], while holographic d-wave models were built up by introducing a charged 
spin two field propagating in the bulk [13, 14]. 

On the other hand, some novel properties of QCD matter in a strong magnetic field 
have been revealed recently, see ref. [15] for a review. In particular, the QCD vacuum 
may undergo a phase transition to an exotic phase with charged p-meson condensed in a 
sufficiently strong magnetic field, which is a kind of anisotropic superconducting phase [16, 
17, 18]. Similar phenomenon has been also studied by holography based on the Sakai- 
Sugimoto model [19, 20, 21]. In particular, the holographic p-meson studied in the confined 
phase at zero temperature [19] suggests that the effective mass of the p-meson would become 
tachyonic in a sufficiently strong magnetic field, which gives rise to the condensate of the 


p-meson in a strong magnetic field. In principle, the holographic p-wave model proposed 
in refs. [11, 12] can also describe the condensate of the p-meson in low energy QCD. 

In a recent paper [11], we constructed a holographic model with vector condensate in 
the Einstein-Maxwell-complex vector field theory in (3 + 1) dimensional spacetime. In 
this model there exists a non-minimal coupling between the complex vector field and the 
Maxwell field. Working in the probe limit, we found that there is a critical temperature 
below which the charged vector condenses via a second order phase transition. The DC 
conductivity becomes infinite and the AC conductivity develops a gap in the condensed 
phase. In addition, it was found that the background magnetic field can induce the conden- 
sate of the vector field even in the case without chemical potential/charge density. In the 
case with non-vanishing charge density, the transition temperature raises with the applied 
magnetic field, and the condensate of the charged vector operator forms a vortex lattice 
structure in the spatial directions perpendicular to the applied magnetic field. Thus, this 
model is also relevant to the p-meson condensate in an applied magnetic field. In a follow 
up paper [12] we studied the model in some detail by taking into account the back reaction 
of matter fields on a black hole background geometry. Depending on the mass and charge 
of the vector field, we found that not only second order, but also first order and zeroth 
order phase transitions exist in this model. In particular the so-called “retrograde conden- 
sation” can also appear, there the hairy black hole solution exists only for the temperature 
above a critical value, its free energy is much larger than the one of the black hole without 
the vector hair. For details please see ref. [12]. 

In the previous studies, we limited ourselves to the black hole background, which corre- 
sponds to a dual field theory in the deconfined phase in a finite temperature, according to 
the gauge/gravity duality. It is worth generalizing our study to the AdS soliton background 
case dual to a field theory in the confined phase. Such a study can not only shed lights on 
the p-wave insulator/superconductor transition at zero temperature, but also is relevant 
to the condensate of p-meson appearing in a confined phase. In this model the background 
magnetic field is from the U(1) gauge field, so we shall call it Abelian magnetic field for 
convenience. The black hole instability induced by a background non-Abelian magnetic 
field coming from a SU(2) gauge field has been reported in refs. [22, 23] where a vortex 
lattice structure forms due to the condensation of some non-Abelian current operators. To 
distinguish our model with the p-wave model with SU(2) gauge field [9], we may call our 
model Abelian p-wave model, while the latter non-Abelian one. To compare two models, 
we will also investigate the AdS soliton instability induced by non-Abelian magnetic field 
in the non-Abelian p-wave model. 

To immerse the strongly coupled system into an external magnetic field background, 
we turn on a uniform magnetic field B in the bulk. For the Abelian p-wave model, we 
find that the increase of the Abelian magnetic field will induce the instability of the AdS 
soliton background. It is well known that the increase of the chemical potential can induce 
the instability of the AdS soliton background, which mimics the insulator/superconductor 
transition in many holographic s-wave [24, 25, 26, 27| and non-Abelian p-wave [28, 29] 
models. The magnetic field here plays the same role as the chemical potential and a family 
of condensate can be induced by the applied magnetic field in the dual strongly coupled 
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system. Although there is a tower of “droplet” solutions in the sense that they are localized 
in a finite region for finite magnetic field [30], it is quite different from the s-wave soliton 
case [31], where the presence of the magnetic field causes the phase transition much more 
difficult. We then study the non-Abelian p-wave model in the AdS soliton background with 
constant magnetic field. Similarly, the non-Abelian magnetic field can also induce the AdS 
soliton instability. In particular, if we take the mass of the charged vector field to be zero 
and the non-minimal coupling constant between the charged vector field and the Maxwell 
gauge field to be one, we find that in the setup of the present paper, the reduced equations 
of motion of the Abelian p-wave model are exactly the same as those of the non-Abelian 
model. In this sense, the Einstein-Maxwell-complex vector field model can be viewed as a 
generalization of the SU(2) p-wave model. We also manage to construct the vortex lattice 
solution in the condensed phase near the transition point, which is very reminiscent of the 
Abrikosov lattice in the common type-II superconductors. 

The plan of this paper is as follows. In section 2, we first introduce the holographic 
model with a complex vector field charged under a U(1) gauge field in the AdS soliton 
background. We turn on a uniform magnetic field and give details about how to recover 
the Landau levels of the vector field in this holographic model. Then we study the phase 
diagram. In section 3, we investigate the magnetic field effect for the non-Abelian p-wave 
model in the AdS soliton background. We construct the general vortex lattice solutions 
in section 4 for both two models. The conclusion and further discussions are included in 
section 5. 


2 Ejinstein-Maxwell-complex vector model 


In this section, we study the AdS soliton instabilities induced by an applied magnetic field 
in the Einstein-Maxwell-complex vector field theory. Let us introduce a charged vector field 
into a (4+ 1) dimensional Einstein-Maxwell theory with a negative cosmological constant. 
The full action reads 


S= 55 | Pav + F- 


1 V 1 V 
qi — sejp” =m pho" + ig7ypupyE™), Q) 
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where K? = 81G is associated to the gravitational constant in the bulk and L is the AdS 
radius which is chosen to be unity. m and q are related to the mass and charge of the 
vector field p,, respectively. F» = V Av — V Ap is the strength of U(1) field A,. The 
tensor p,, is defined by py» = Dypy — Drp, with the covariant derivative D, = Vp — iqApn. 
The last term in (1) describes the non-minimal coupling between the charged vector field 
Py and U(1) gauge field A, characterizing the magnetic moment of the vector field p,,. The 
coupling constant y in some sense can be thought of as an effective gyromagnetic ratio of 
the charged vector field. 

The equation of motion for the charged vector field p, is given by 


D ig = m Pye Fy, = 0, (2) 


while the equation of motion for gauge field A,, is 


VF, = igl ey — P Pup) + igy” (pvp, — pl py). (3) 


By taking the limit q — œ with qp„ and qA, fixed, the back reaction of the matter sources 
to the background spacetime can be ignored. This is the probe limit we will adopt in this 
paper. Since we want to study the phase transition in the confined phase, the background 
geometry is taken to be a (4 + 1)-dimensional AdS soliton 
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ds? = a + r?(—dé? + dr? + dy’) + f(r)dy’, (4) 
where f(r) = r?(1 — 2) and rọ is the tip of the soliton. This is a self-consistent solution 
of action (1) by neglecting the matter sources. To avoid the potential conical singularity 
at r = ro, the spatial direction 7 must have a period n ~ 7+ mT/ro. The AdS soliton just 
looks like a cigar with the asymptotical geometry R!? x St near the AdS boundary. Since 
the spacetime exists only for r > ro, the dual field theory is in a confined phase and has 
a mass gap, Eg ~ ro. The AdS soliton spacetime can describe an insulator in condensed 
matter physics [24], due to the presence of the mass gap. 


2.1 Turning on a constant magnetic field 


We now turn on a magnetic field to study how the applied magnetic field influences the 


system. Similar to the black hole case studied in ref. [11], a consistent ansatz turns out to 
be 


pdx” = [ep,(r, £)” + O(e?)Jdx + [epy(r, zee” + O(c?) dy, P 
Adz” = [ġ(r) + O(®]dt + [Bx + O(e?)]dy, (5) 
where pz(r, £), py(r, x) and ¢(r) are all real functions, p is a real constant and the constant 0 
is the phase difference between the x and y components of the vector field p,,. The constant 
Abelian magnetic field B is perpendicular to the x — y plane. The parameter € is a small 
quantity characterizing the deviation from the critical point at which the condensate begins 
to appear. 
Substituting the ansatz (5) into equation (3), the equation of motion for ¢ can be read 
off at the zeroth order 
Pal 


+e) =0. (6) 


The asymptotic value of ¢ gives the chemical potential u = A;(r — oo) of the dual field 
theory. Imposing the Neumann-like boundary condition at the tip of the AdS soliton 
background [24], we find a unique solution 


lr) = pb. (7) 


o'(r) + ( 


The equations of motion for pẹ and py can be deduced from (2) at linear order O(e). We 
can further make a variable separation as p,(r,x) = y,(r)X (x) and p,(r, x) = py(r)Y (2). 
To satisfy the equations of motion, @ can only be 64 = 5 + 2nm or 0- = —5 + 2nr with n 
an arbitrary integer. We then obtain the equations of motion for y,(r), y,(r), X(x) and 
Via) as 


y,X + (qBx — p)p,Y =0, (8) 
pX + (qBx — p)p,Y =0, (9) 
ote (Lp dyyt 4 (Ce yp. 22 fe (qBe—p) aa 2 (gBe—p)"] = 0, (10) 
ay ae eee) E a | X Dz X í 
T PR m? G X p: X pa 
EN a ee le EN E E E B20 (i 
ey her ey a 7? yt rafly (qBx yo, (1+7) rae , (11) 


where the prime and dot denote the derivatives with respect to r and x, respectively. Here 
and below the upper signs correspond to the 0, case and the lower to the 0_ case. To 
satisfy (8) and (9), one should impose the constraint 


Py =Chn, X+c(qBu—p)Y =0, (12) 


with c a real constant. We can see that only two of the four functions are independent. 
Substituting (12) into the remaining equations, we can find the following three equations 


fio Qu -E m 

" J a. UP Te He — 

— X +c(1+y)qBY + (qBz — p)’ X = EX, (14) 
=F ey + (qBz — p)?Y = EY, (15) 


where € is a constant. It is clear that (14) and (15) are the same as corresponding ones 
in the black hole case discussed in ref. [11]. For the sake of brevity, we therefore just list 
the main results in this paper. More details can be found in ref. [11]. By defining a new 
function as 


ylz) = X(x) — Y(2), (16) 
one gets its equation from (14) and (15) as 
Y(z) + [E F c(1 + 7)qB — (B7 — p)"}b(2) = 0, (17) 


where we have imposed the condition c? = 1 [11]. This eigenvalue equation can be solved 
exactly with the eigenvalue function 


Yalt) = Nye 28\@- a5)" Hy, (4/|q BI (a — B (18) 


and the corresponding eigenvalue 
En = (2n + 1)|4B| + (1 + 7)4B, (19) 


where H,, is Hermite function, N, is a normalization constant and n is a non-negative 
integer labeling each “energy level”. 

The solutions of Y, and p} corresponding to the eigenvalue €,,, denoted by Yzn and Pyn, 
can be obtained by solving the equation of motion (13) with €, given in (19). So far, we 
have recovered the Landau levels. As we arrive at the transition point, we should encounter 
a marginally stable mode at the critical point. Those solutions obtained from (13) just 
correspond to the marginally stable states. 

One can see from (13) that the effective mass of p, is given by 


2 En- Gi a (2n 1) eB sacl +yaB -Pp 
2E mt O LOUE, 
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= - (20) 
which is clearly corrected by the magnetic field B. The appearance of magnetic field can 
increase or decrease the effective mass, thus will hinder or enhance the transition from the 
normal phase to the condensed phase. In what follows, we consider the solution with the 
lowest Landau level with n = 0, which reads 


& = —|y4B], 
No aB (21) 
X¢(a;p) = Se Oa” = -Yy (eip). 
The general falloff of pẹ near the boundary r —> oo behaves as 
Px Pr+ 
£ = ee 22 
a ae pk (22) 


with A4 =1+V71+ m?. According to the AdS/CFT dictionary, up to a normalization, 
px and pr, timing the spatial dependent part are regarded as the source and the expec- 
tation value of the x component of the dual vector operator denoted as J” in the boundary 
theory, respectively. Similar statement holds for y, corresponding to the y component of 
the dual vector operator J¥. Since we want the U(1) symmetry to be broken spontaneously, 
we turn off the source term, namely take p,_ = 0. 


2.2 Phase diagram 


We are interested in how the applied magnetic field can influence the critical point from 
the normal phase to the condensed phase. As we can see from (20), the effective mass of 


'The m? has a lower bound as m? = —1 with A, = A_ = 1. In such a case, there is a logarithmic 
term in the asymptotical expansion. We treat such a term as the source which is set to be vanishing to 
avoid the instability induced by this term. 


the charged vector field p, in the lowest Landau level n = 0 depends on the magnetic field 
B and the non-minimal coupling parameter y as 
B 25,2 
Mee =m — lyaB| + g'u" (23) 


fF 
e r2 


It is clear that the increase of the magnetic field B decreases the effective mass and thus 
tends to induce the phase transition. The magnetic field plays the same role as the chemical 
potential ju. 
We introduce a new coordinate z = ro/r, then the equation of motion (13) can be 
rewritten as °? 
14324 , m? A 


plz) — za wre) = laa E Tael =0, (24) 


. We further define a new function for numerical convenience 


2,2 
B 
where A = vu hah 
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Figure 1: Three marginally stable modes of the charged vector field. The three curves 
from top to down correspond to Ag œ 7.7562 (red), Ay œ 27.992 (blue) and A, œ 59.601 
(purple), respectively. We here take m? = 5/4. We have normalized the value of Y, at the 
tip to be one. 


palz) = (°F) (25) 


then the equation of motion for F(z) is given by 


7 (a — 2A_)F'(2) — n = A-(A- +2) — E =. (26) 


F"(z) 


2 
2For brevity, we also use prime to denote the derivative with respect to z. But the meaning of the 
derivative in the text is clear and will not be confused. 


To solve such a second order equation, we should impose suitable boundary conditions. 
Remember that the leading term of F near the boundary z — 0 gives the source term 
which is set to be zero, so one has the condition F(z = 0) = 0 at the boundary. We 
impose the regularity condition at the tip z = 1. Thanks to the linearity of (26), we take 
F(z = 1) = 1 in our numerical calculation without loss of generality. It is clear that (26) 
depends on two parameters m? and A. For a given m7, only for certain values of A can the 
boundary conditions be fulfilled. 

Solving (26) by shooting method, we present the three lowest-lying marginally stable 
modes in figure 1 for m? = 5/4. The three modes are in the sequence Ag < Ay < As. The 
mode with Ag ~ 7.7562 (the red curve) is regarded as a mode with node n = 0 since the 
curve has no intersecting points with horizontal axis except at the origin z = 0. Following 
such terminology, the blue curve corresponding to A; œ 27.992 is considered as a mode 
with node n = 1 and the purple one corresponding to Ag œ 59.601 as a mode with node 
n = 2. Due to the radial oscillations in z-direction of y,(z), the latter two modes are 
therefore thought to be less stable than the first one. Thus, the lowest value Ap just gives 
the critical value above which the AdS soliton is unstable to developing a vector “hair”. 
We plot the value of Ag for different squared mass of the vector field in figure 2, from which 
one can see that Ap increases as we increase the squared mass. 


Figure 2: The critical value of Ag with respect to m? of the vector field. The black points 
are obtained by solving the equation (26) numerically by use of shooting method. 


It should be stressed here that it is the combination Ag = Petan] that determines the 
phase transition point of the system. More specifically, turning off the chemical potential 
u, we can see that only the magnetic field itself can also induce the instability. This 
result is totally different from the case in the holographic s-wave superconductor/insulator 
model [31], where the presence of the magnetic field makes the normal phase more stable 
rather than unstable. The result here is reminiscent of the QCD vacuum instability induced 
by strong magnetic field to spontaneously developing the p-meson condensate. Just as 


pointed out in the finite temperature case [11], this is due to the non-minimal coupling 
of the vector field p” to the gauge field A, in the last term in (1). The reader would be 
familiar with such a term used to describe the coupling of magnetic moment for charged 
vector particles to a background magnetic field in the literature, see refs. [32, 33] for 
example. 


3 Ejinstein- Yang-Mills model 


The black hole instability induced by an applied non-Abelian magnetic field in the SU(2) 
p-wave model was studied in refs. [23, 34]. In this section, we will study the AdS soliton 
instability induced by the non-Abelian magnetic field. 

The five-dimensional Einstein-SU(2) Yang-Mills theory with a negative cosmological 
constant can be written as [9] 


1 a EE 
S= | &xv=aggR+ 7) 49 res a J, (27) 


where ĝ is the Yang-Mills coupling constant. The SU(2) Yang-Mills field strength is given 
by 
Fe = OA, = DA, + ora A (28) 


preys 
123 _ 


where e% is the totally antisymmetric tensor with e! = +1. The gauge field is given by 


A= Ajr*dz". The equation of motion for Aj, reads 
Vo TA Fer =o (29) 


We work in the probe limit and treat marginally stable perturbations as probe into the 
AdS soliton (4). 


3.1 Adding a constant magnetic field 


We turn on a background non-Abelian magnetic field B and also allow for a non-Abelian 
chemical potential in the dual system. Both of them come from the 7° component. of the 
Lie algebra. We take the following ansatz 


Alda" = [ea}(r, x,y) + O(€)]dx + [ea] (r, x,y) + O(€)]dy, 
Ade!" = [ea?(r, x,y) + O(€)]dx + [ea?(r, x,y) + O(€)]dy, (30) 
Ade" = [6(r) + O(2)|dt + [Be + O(2)]dy, 


with e the small deviation parameter. Substituting above ansatz into (29), we obtain the 
equation of motion for ¢(r) at leading order, which is just the same as (6) with the solution 
given in (7). At linear order O(c), the equations of motion turn out to be 


sVs + (Oy —iBx)V, = 0, (31) 


3u,v = (t,r,2,y,n) denote the indices of spacetime and a,b,c = (1,2,3) are the indices of the SU(2) 
group generators T° = 0% /2i, where o% are Pauli matrices. 
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ðzð-Vz + (Oy — iBx)d,V, = 0, (32) 


/ 
BVHE EOV So [(0fi2B10,-+ 4° —B*2* Vert (-0,04-+4B20, iB Vy = 0, (33) 
2 F 1 1 : * 2 2 
3V; + (F + =)3-Vy + pap ley +iBxô, + i2B)V; (0, + )Vy| = 0, (34) 


where V, = a} — ia? and V} = a, — ia’. 
In order to solve above four equations of motion, we further make a variable separation 


V,(r,2,y) = Galr)X(a)e™, Vy (r,2,y) = Gy (r)Y (zje, (35) 


where Gz, Py, X and Y are all real functions and p is a real constant. We also introduce 
a constant @ describing the phase difference between V, and V,. Further analysis shows 
that 6 can only be 6, = 5 + 2nm or @_=-= 5 + 2nn with n an eee integer. Then we 
obtain the following equations of motion 


GX + (Ba — p)6,¥ =0, (36) 
G,X + (Ba — p)@Y =0, (37) 
1 2 p- Te tare 
<" Í dee H ~ Pu y y Y y 2 
ğ (F )e ap? paper are X Ge ( p)] (38) 
sic af. Leen Ê. gY XG, X bp 
(pE dG le xr +2B = 0, 39 
Py (5 ey T apy ( FG, To, ( ) 


where the prime and dot denote the derivatives with respect to r and x, respectively. The 
same as in the previous section, the upper signs correspond to the 6, case and the lower 
to the 6_ case. Comparing the above four equations of motion with equations (8)-(11), 
we can find that the two sets of equations of motion are the same if we choose parameters 
q = 1,m? = 0 and y = 1 in (8)-(11). Thus, the results obtained in the previous section 
can be used here by simply restricting the model parameters to q = 1,m? = 0 and y = 1. 
In this sense, the complex vector field model can be thought of as a generalization of the 
SU(2) p-wave model. 


3.2 Phase diagram 


In the presence of constant magnetic field and chemical potential, the field y, * acquires 
an effective mass in the lowest Landau level n = 0 as 


2 Eo — E |B| +p’ 
me = — z — a eT 


(40) 


r 


4Since the non-Abelian p-wave model in our ansatz is the same as the Abelian p-wave model with 
parameters q = 1,m? = 0 and y = 1, in what follows, we shall omit the tilde in all functions in this 
non-Abelian model for clarity. 
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It is clear that the increase of the magnetic field B decreases the effective mass and thus 
tends to induce the instability. For given chemical potential u, as we increase the mag- 
netic field B to the critical value Agr§ — u? with Ao = (u? + |B|)/rő ~ 5.1313, the AdS 
soliton background will become unstable to developing non-trivial vector “hair”. It is clear 
here that even without the chemical potential, the magnetic field will induce the vector 
condensate and the AdS soliton instability. In addition, note that the critical value Ag is 
comparable to the critical magnetic field found in the black hole background [34]. 


1.0f 
0.8 i condensed phase | 


0.6} 
/Mcl | 
0.47 normal phase 


0.2} 


Ge ee ee a, 
0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 3: The value of chemical potential from the normal phase to the condensed phase as 
a function of magnetic field. It corresponds to the lowest Landau level with Ap œ 5.1313. 
He = v Aoro is the critical chemical potential in the case without the applied magnetic 
field. 


The (u, B) phase diagram for the lowest Landau level with Ao ~ 5.1313 is drawn in 
figure 3. In order to determine which side of the phase boundary is the condensed phase, 
we take a look at the equation (40) which suggests that the magnetic field decreases the 
effective mass. As one increases the strength of magnetic field at a fixed chemical potential, 
the normal phase will become unstable for sufficiently large magnetic field. Therefore, the 
region in the upper right in figure 3 corresponds to the condensed phase. In addition, let 
us mention here that the (u, B) phase diagram for a fixed m? for the Abelian p-wave model 
discussed in the previous section should be extremely similar to figure 3. 


4 Vortex lattice solution 


Following ref. [11], let us now construct the vortex lattice solution for the Abelian and the 
non-Abelian p-wave models. We consider the case with the lowest Landau level n = 0. 
Since the eigenvalue Eo is independent of p, it is easy to find that a linear superposition of 
the solutions e” pro(r) Xë (x; p) and e”¥y,o(r) Yo (x; p) with different p is also a solution 
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of two models at linear order O(€). Therefore, the following two functions 


+00 
W,(r,2,Y) = par) S cee” X5 (x; pe), 
L=—00 
W,(r, x,y) = ce™*pao(r) A cee¥Y (a; pe), 
lL=—0o 
. . . o iE B lqBle . 
are solutions of the full equations of motion, where cg =e “îi and pẹ = with 


a, and az arbitrary constants. According to ref. [11], we can construct the vortex lattice 


solution 


AW (r, x,y) = W,(r, £, y) — ce =W, (r, £, y) = prolr 5 cee” aby (x; pe), (42) 


L=—00o 


where polz; pe) = Noe ~210Bl(e— a5)” 
The solution AW (r, x, y) exhibito a pseudo-periodicity 


AW ga =AW mna, 


vlaB[ 


43 
AW + a + eg) TVHD AWG 259) = 
and has a zero at 1 i 
Xin. = (M+ 5)vi + (ñ + 5)vo, (44) 
where two vectors vı = aia and vg = PT + nei! [35]. Mm and ñ are two 


integers. Remember that the corresponding results for the SU(2) model can be directly 
obtained by choosing the parameters q = 1,7 = 1 and m? = 0. 

According to the AdS/CFT dictionary, the coefficient of the sub-leading term near the 
boundary r + oo of AW (r, x, y) just gives the vacuum expectation value of the dual vector 
operator. For the Abelian p-wave model, the expectation value of the operator Je dual to 
pz is given by the coefficient of 1/r4+ at boundary r > oo, the combination J} = (Jei J») 
dose exhibit the vortex lattice structure which has the cores located at Xma. This result 
is the same as that we reported in the black hole background [11]. In particular, to obtain 
the triangular lattice we can choose the following parameters 


2 2 
“es 2v7 ae, (45) 


V3” v3 
We notice that it is the linear combinations J} which exhibit the vortex lattice solution. 
In particular, for q > 0, it is J_ that corresponds to the lowest Landau level, while for 
q < 0, it is J}. The SU(2) p-wave model corresponds to the case with q = 1, so its vortex 
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lattice solution is related to J_. This further provides the evidence for the correctness of 
choosing c? = 1 in section 2.1. 

As a concrete example, the vortex lattice is presented in figure 4 in the x — y plane 
perpendicular to the magnetic field for the triangular structure. For other values of a; and 
a2, one can obtain different vortex lattice configurations. To determine the true ground 
state, one must go beyond the linear analysis to find the values of a; and a2 which minimize 
the free energy of all possible lattice structures. What’s more, the overall normalization of 
AW (r, x,y) can only be fixed at higher order. The calculation of higher order contributions 
can be done straightforwardly as in refs. [22, 35]. Nevertheless, it is much more involved 
and is not very relevant to our purpose of this paper. Thus, we shall leave it for our further 
study. 
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2 4 
Figure 4: The contour plot of the vortex lattice structure for the triangular lattice in 
x — y plane. Darker shading indicates smaller values of the norm of the condensate. In 
particular, the condensate vanishes at the core of each vortex. 
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5 Conclusion and discussion 


As a continuation of our recent studies on a holographic p-wave superconductor model in 
the Einstein-Maxwell-complex vector field model proposed in ref. [11], in this paper we 
generalized to the case with vector condensate of dual conformal field theory in confined 
phase, i.e., in an AdS soliton background (see the metric (4)). This study is relevant 
to the p-wave superconductor(superfluid) /insulator phase transition as well as charged p- 
meson condensation in confined QCD vacuum in strong magnetic field. We focused on the 
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behavior of the model in the presence of a constant magnetic field B. At linear order O(e), 
the spatial part of vector field p,, can be solved analytically, which gives the Landau level 
of the vector field in a magnetic field background. It was found that the effective mass of 
the vector field is reduced by turning on magnetic field in the lowest Landau level case. 
More precisely, it is the combination A = (q?u? + |yqB|)/ré that determines the instability 
of our system, here u is the chemical potential of dual field theory, while rọ is related to the 
mass gap of the field theory. When we increase the value of A to the critical one, the AdS 
soliton background becomes unstable to developing non-trivial vector “hair”. Interestingly, 
it can be seen from A that the magnetic field itself can induce the instability even without 
chemical potential u. It is clear that the essential point for this is due to the non-minimal 
coupling between the vector field and the background magnetic field. 

To have further understanding, we also study the Einstein-SU(2) Yang-Mills theory 
in the AdS soliton background with a non-Abelian magnetic field. Interestingly, at least 
under our ansatz in the present paper, we found that the reduced equations of motion for 
the non-Abelian model are the same as those of the complex vector field model in the case 
with parameters m? = 0 and y = 1. ° Therefore, the presence of a non-Abelian magnetic 
field can also make the AdS soliton background unstable. In some sense, the Einstein- 
Maxwell-complex vector field model is a generalization of the SU(2) model to a general 
mass square m? and gyromagnetic ratio y of the vector field. 

The phase diagram in terms of B and pu is drawn in figure 3. The critical chemical 
potential above which a condensed phase appears decreases with the applied magnetic 
field, which means the presence of magnetic field can enhance rather than suppress the 
phase transition. This result is reminiscent of the condensation of W-meson [36] and p- 
meson [16, 17] in a strong magnetic field. Of course, this result is opposite to the well 
known behavior reported in ordinary superconductors as well as holographic s-wave su- 
perconductor/insulator model [31], there the presence of an applied magnetic field will 
suppress the superconducting phase transition. Remember that the action (1) describes a 
p-wave model, and that for the p-wave superconductors, theoretical studies [37, 38, 39, 40] 
as well as experimental evidences [41, 42] suggest the possibility that the magnetic field 
can induce superconductivity in some cases, in particular, for the ferromagnetic supercon- 
ductors. Therefore our results are understandable. 

In the present study, we limited ourselves to the probe approximation by neglecting 
the back reaction of matter fields on the bulk AdS soliton geometry, which may work 
only near the critical point with continuous phase transition. Although it does be able to 
reveal some significant properties, something would be lost in this approximation, especially 
for the whole phase structure, see ref. [43] as an example. Indeed, in the black hole 
background case, going beyond the probe approximation, we found a rich phase structure in 
this Einstein-Maxwell-complex vector field model [12]. Depending on mass square m? and 
charge q, the phase transition can be zeroth order, first order or second order. In addition, 
there also exists the so-called “retrograde condensation”. The behavior of entanglement 


©The condition q = 1 is not important, since the charge q of the vector field p, can be scaled to be one 
in the probe limit. 
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entropy in this model has also been studied in ref. [44]. In some sense, although the 
Einstein-Maxwell-complex vector model looks to be a generalization of the SU(2) p-wave 
model, it can have much rich phase structure and distinguished properties compared to the 
latter. Therefore, it is interesting to go beyond the probe limit in the soliton background. 
Our ongoing work shows that this model has much more interesting behaviors in the soliton 
background. 

The Einstein-Maxwell-complex vector field model described by the action (1) is just a 
phenomenological model. At the moment it is not sure whether this model can be embedded 
into some effective theory of string/M theory. If the answer is yes, it is certainly of great 
interest to understand those interesting properties of the model from the dual field theory 
side. This must be quite helpful to understand some high temperature superconductors 
and QCD vacuum in strong magnetic field. 
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